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ABSTRACT 



Context, d Velorum is a nearby (~24pc) triple stellar system, containing a close, eclipsing binary (Aa, Ab) discovered in 2000. 
Aims. Multiple systems provide an opportunity to determine the set of fundamental parameters (mass, luminosity, size and chem- 
ical composition) of coeval stars. These parameters can be obtained with particular precision in the case of eclipsing binaries; for 
6 Velorum's components (Aa, Ab) this potential has however not yet been exploited. 

Methods. We have analyzed interferometric observations of the close binary (Aa, Ab), obtained with the VINCI instrument and two 
VLTI siderostats. The measurements, which resolve the two components for the first time, are fitted onto the simple model of two 
uniformly bright, spherical stars. 

Results. The observations suggest that Aa and Ab have larger diameters than expected if they were on the main sequence, and that 
they are, thus, in a later evolutionary state. 

Key words, multiple stars - interferometry - stellar evolutionary state 



1. Introduction 

5 Velorum (HD 74956) is one of the fifty brightest stars on the 
sky, with a visual magnitude of m v = 1 .96 mag (Johnson et al. 
1966). It is a multiple stellar system (e.g. Worley & Douglass 
1 19971 1, but in spite of its brightness and near distance, n = 
(40.90 ±0.38) mas (Perryman et al. |1997l ), the issue of its compo- 
sition remains unresolved. As early as 1847, Herschel published 
his detection of two faint visual companions, 6 Vel C and D, at 
distance 69" apart from 5 Vel A. Another companion, 5 Vel B 
- at the time separated by ~ 3" from 6 Vel A - was later dis- 
covered by Innes J1895I ). The separation 07736 + 07014 between 
the components A and B appeared surprising when measured by 
Hipparcos, but it was well explained later, in terms of the orbit 
computation of Argyle et al. (2002) that showed a highly ellipti- 
cal orbit of component B with period P = 142 yr. In 1979, pre- 
liminary results from speckle interferometry suggested an even 
further component of the system (Tango et al. |1979b . This appar- 
ent companion was found at a separation of ~ 076 and was taken 
to be a further component, because for star B the separation at 
the time was believed to be ~ 3". By now, however, it seems very 
likely that the speckle observations resolved 5 Vel B; while there 
is still an unexplained disagreement for the position angle, the 
measured small separation does fit well with the orbital solution 
found by Argyle et al. 020021 ) . As noted in earlier publications, 
(Hoffleit et al. |1991| Otero at al. 120001 ), the two stars that are cur- 
rently termed 5 Vel C and D were taken to be associated with the 
pair AB because of seemingly similar proper motion. However, 
we have not found the source of the proper motion measurement 
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of C and D. Finally, the most luminous component, A, was re- 
cently recognized to be a close eclipsing binary with a period 
T = 45.15 days (Otero et al. 2000). 6 Vel has since been classi- 
fied as a quintuple stellar system. 

This investigation is focussed on the bright eclipsing binary, 
6 Vel A, but we also argue that 6 Vel C and D are not physically 
associated with £ VelA,B. While this makes 6 Vel a triple sys- 
tem, it takes away little from its challenging potential to obtain 
important information on stellar evolution. As the inclination, i, 
of its orbital plane is constrained to be close to 90°, an eclipsing 
binary system provides one of the best means to obtain, in terms 
of the Kepler laws of motion, fundamental stellar parameters. 

In this research note we present the first interferometric ob- 
servations of the eclipsing binary 6 Vel A, obtained with ESO's 
Very Large Telescope Interferometer (VLTI) and its "commis- 
sioning instrument" VINCI. The measurements resolve, for the 
first time, this binary system. They are here analyzed with non- 
linear least-square fitting methods. We combine our interfero- 
metric results with existing photometric and spectroscopic ob- 
servations, estimate some orbital parameters of the S Vel A bi- 
nary system, and discuss, based on the results, the stellar prop- 
erties of the individual components. 



2. Characteristics of 6 Vel A derived from previous 
measurements 

2.1. Orbit orientation and eccentricity 

In this subsection, a priori estimates of two orbital parameters 
of the 6 Vel (Aa-Ab) system are derived from the time interval 
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Fig. 1. The angle, u>, at primary eclipse and the eccentricity, e, 
are constrained by the fractional durations between the eclipses, 
Tf - 0.43 ± 0.05, and the fractional durations of the eclipses, 
p f = 1.78 ±0.19. 



between the eclipses and their durations, the eccentricity, e, and 
the angle, w, between the semi-major axis and the line of sight. 
The angle a> is similar to, but must not be confused with the more 
generally used parameter longitude of periastron. 

As reported by Otero et al. ([20001 and Otero (|2006|>, the frac- 
tional orbital period from the primary to the secondary eclipse 
equals r/ = 0.43 ± 0.05. The secondary eclipse has been ob- 
served by the Galileo satellite in 1989, and its duration and depth 
have been fairly precisely established to be 0.91 ± 0.01 days and 
Am// = 0.32 ± 0.02 (Otero |2006] l. The same spacecraft observed 
the primary eclipse several years later, although its measure- 
ments had become less accurate then. The approximate dura- 
tion and depth of the primary eclipse are 0.51 ± 0.05 days and 
Am/ = 0.51 ± 0.05 (Otero 2006). The ratio of durations thus 
amounts to p f = 1.78 ± 0.19. 

The relative motion of the two stars 6 Vel Aa and Ab is taken 
to be independent of external forces, and the vector, s, from Ab 
to Aa traces an elliptical orbit around Ab as a focal point. As 
the photometric light curve indicates a total eclipse for 6 Vel A, 
the inclination of the orbit needs to be close to 90°. To simplify 
the equations, we assume i = 90°. Ab is taken to be the star 
with the higher surface brightness. During the primary eclipse, 
which is deeper, Ab is thus eclipsed by Aa. The angle of s, 
also called the true anomaly, is zero at periastron and increases 
to n as the star moves towards apastron. The distance between 
the stars depends on according to the relation: 



s(0) = a(l -e 2 )/(l + ecos(0)) 



(1) 



where a denotes the semi-major axis. In line with Kepler's sec- 
ond law the vector s covers equal areas per unit time. The frac- 
tional orbital period to reach angle is, accordingly: 



(2) 
(3) 



Tiff) = \ f S 2 (0')dff' 

A Jo 

= [2arctan(/i tan(0/2) - f 2 sin(0)/(l + e cos(0))]/27r 

where A = na 2 - e 2 ) equals the area of the ellipse. f\ = 

V((l - e)/(l + e)) and f 2 = e ^j(jTe ^ ). 

During the primary eclipse, when the star with the lower sur- 
face brightness, Aa, covers Ab, the vector s is directed towards 
Earth, and 6 equals a>. During the secondary eclipse 6 equals 



a) + n. Thus: 

Tf = 0.43 ±0.05 

r 

d0/(l +ecos(0)) 2 



-f 



(4) 
(5) 



= [arctan(/i tan((w + n)/2)) - arctan(/i tan(w/2)) 

+ f 2 sm(cj)/(l-e 2 cos(aj) 2 )]/n (6) 

which determines co for any given eccentricity, e (Fig.Q]). The 
orbital velocity decreases as goes from zero to n, i.e. from the 
periastron to the apastron. In the subsequent interval, n to 2n 
(or -7T to 0), it increases again. If the line of sight contained the 
orbital major axis, i.e. a> = or n, the fractional duration be- 
tween eclipses t would equal 0.5. Note that for such values of u, 
Eq.|6]is not defined, yet r tends towards 0.5 when u approaches 
0, or n. If the line of sight contained the orbital minor axis, i.e. 
cj - n/2 or -n/2, the maximum and minimum values of r would 
be reached. Values of t less than 0.5 are thus associated with 
negative u> values. Since the fractional orbital period from the 
primary to the secondary eclipse is 0.43, the angle u must lie 
between -n and 0. As Fig. Q] shows, the eccentricity needs to be 
larger than m 0.03. 

On the other hand, a> can be further constrained through the 
ratio of the eclipse durations as follows. The eclipse durations 
are inversely proportional to the product rd0/dt of radius and 
angular velocities during the eclipses. They are thus proportional 
to s(6), and their ratio is: 



p(u>) = (1 - ecos(<y))/(l + ecos(w)) 



(7) 



Given pf = 1.78 ± 0.19, this leads to a second relation between 
e and u>. As illustrated by Fig.Q] simultaneous agreement with 
both observed values Tf and pf is reached only if e e [0.23- 
0.37] and w e -[0.1 -0.7]rad. 

2.2. Semi-major axis and stellar parameters 

Orbital motion in the triple system 5 Vel(Aa+Ab+B) has re- 
cently been substantiated and analyzed by Argyle et al. (|2002). 
From position measurements taken over a period of roughly 
100 years the authors inferred a P = 142 yr orbit for component 
B, and deduced a total dynamical mass M(Aa)+M(Ab)+M(B) = 
5.7*}qo M e . Photometric and spectroscopic measurements of the 
individual components being few and partly inconclusive, indi- 
vidual mass estimates are still difficult. 

Hipparcos measured an apparent magnitude of Hp = 1.991 
for 5 Vel A and H P - 5.570 for 6 Vel B. With the transformations 
given by Harmanec (U998I I, the approximate Johnson V mag- 
nitudes are my = 1.99 and my = 5.5 for 6 Vel A and d> VelB, 
respectively. But the colours of the individual 6 Vel components 
being unknown, it needs to be noted that the uncertainty of my 
can be as high as ~ 0.07mag. Since 5 Vel is close (d = 24.45 pc 
according to Hipparcos) no interstellar reddening towards the 
source needs to be assumed, and hence, the absolute magnitudes 
are My ~ 0.05 for 6 Vel A and My ~ 3 .6 for 6 Vel B . 

Several authors have analyzed spectra of 6 Vel A (e.g. Wright 
120031 Alekseeva [19971 Levato [19721 Gray & Garrison \TWf\ . 
Many of their measurements have probably included i5VelB, 
but its flux is too low to add a significant contribution. From 
the metal line ratios and Balmer line equivalent widths all au- 
thors deduced either spectral type AO V or Al V. This being most 
likely an average classification of the two stars, Aa and Ab, one 
star should be slightly hotter and the other cooler than a AO/ IV 
star. No signatures of a double-lined spectroscopic binary were 
reported in any of the spectroscopic observations. 
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Based on the spectrophotometric information referred to 
above, and under the assumption that all 6 Vel components are 
on the main sequence, it is suggested that Aa and Ab are of spec- 
tral type between AOV and A5V with masses in the range 2.0- 
3.OM . Furthermore it follows that B is an F-dwarf with mass 
about ~ 1.5M . This agrees reasonably with the total dynamical 
mass derived by Argyle et al. (2002]). 

An a priori estimate of the semi-major axis, a, of the Aa-Ab 
system is next derived from the mass sum of Aa+Ab (5 ± 1 M ) 
and its orbital period (T = 45.150 + 0.001 days), which leads to 
a = (6.4 + 0.5) x 10 10 m=0.43 + 0.04AU. 
If they are main sequence early A stars, Aa and Ab should have 
stellar diameters between 1 .7 - 2.4 D Q . 

Finally, the depths of the eclipses can be used to constrain 
the surface brightness ratio <p of the two eclipsing components, 
6 Vel Aa and Ab, 



1.28 < ■ 



1 _ lQ-Am,/2.5 

I _ io- A "'"/ 2 - 5 



< 1.67 



(8) 



3. VLT Interferometer/VINCI Observations 

3.1. Data description 

During April-May 2003, the ESO Very Large Telescope 
Interferometer (VLTI) was used to observe the eclipsing bi- 
nary 5 Vel( Aa+Ab) in the K-band at four orbital phases with the 
single-mode fiber based instrument VINCI (Glindemann2001 , 
Kervella et al. 2003 ). The observations were performed with two 
siderostats, placed at stations B3 and M0, separated by 155.368 
m. TableQ] summarizes the resulting data. 
Every interferometric observation yields a fringe contrast or 
squared visibility, V 2 , whose variations are not only due to inter- 
ferometric modulation, but also to atmospheric and instrumen- 
tal fluctuations. Accordingly the raw squared visibilities need to 
be calibrated by a reference star. To this purpose the observa- 
tions of 6 Vel were combined with observations of HD 63744, 
a star of spectral type K0III, with an estimated diameter of 
1.63 + 0.03 mas (Borde et al. 120021 . The interferometric mea- 
surements were then analyzed by use of the VINCI data reduc- 
tion pipeline, described in detail in Kervella et al. ( 120031 1. 

Additionally, the calibrated V 2 values need to be corrected 
for the influence of the nearby component (TVel B. The diffrac- 
tion on the sky (through an individual VLTI 0.4 m siderostat) of 
the fiber fundamental mode, which defines the interferometric 
field of view, is equivalent to an Airy disk with a 1'.'38 diameter. 
At the time of the observations Aa+Ab and B were separated 
by ~ (1.0 ± 0.3)". Depending on atmospheric conditions, the in- 
terferograms are, therefore, contaminated by a random and time 
varying fraction of light, i.e. an incoherent signal, from star B. 
The visibilities must, accordingly, be multiplied by a factor: 



V c = V X (1 + / B //Aa + Ab) = (1-05 ± 0.05) x V 



(9) 



/b and /Aa+Ab are the intensities collected by the interferome- 
ter from £Vel B and <SVel (Aa+Ab). Ig/lAa+Ab lies between (no 
light from B) and 10~ Am/2 5 = 0.09 (star B is completely in the 
field of view), where Am ~ 2.6 equals the K-band magnitude 
difference between B and Aa+Ab. 



3.2. Comparison to a model 

The 17 visibility measurements, V 2 , were fitted to a model of a 
binary system of two uniformly bright spherical stellar discs, ob- 
served at K-band with a filter of finite bandwith. Five parameters 



Table 1. Details of the VINCI measurements: observing date, 
orbital phase of 6 Vel (Aa+Ab), calibrated squared visibility V 2 
and standard deviation cry 2 , number of accepted scans Ns (out of 
500). The uncertainty on the phase determination equals ±0.002. 



Date 


Julian Date - 2452700 


Phase 




N s 


21 Apr 03 


50.628 


0.937 


57.40+3.60 


383 




50.633 


0.937 


54.20+3.60 


298 




50.639 


0.937 


54.00+3.50 


311 


03 May 03 


62.498 


0.200 


27.54+0.66 


96 




62.502 


0.200 


34.03+0.70 


393 




62.507 


0.201 


43.40+2.20 


260 




62.512 


0.201 


42.20+5.07 


68 




62.542 


0.201 


13.37+0.45 


80 




62.545 


0.201 


8.47+0.58 


122 




62.554 


0.202 


5.06+0.17 


356 




62.562 


0.202 


15.32+0.33 


416 


10 May 03 


69.551 


0.357 


44.30+1.80 


435 




69.556 


0.357 


52.20+2.00 


446 




69.561 


0.357 


56.40+2.10 


455 


1 1 May 03 


70.492 


0.377 


8.30+0.45 


258 




70.506 


0.378 


2.92+0.40 


116 




70.519 


0.378 


1.30+1.40 


45 




Fig. 3. x 1 as a function of the stellar diameters. The three other 
parameters of the model are set equal to: a = (5.7 ± 0.3)10 10 m, 
e = 0.230 ± 0.05, D. = 27 A + 1.2°. 



of the binary model (stellar diameters D a , D/,, position angle of 
the Ascending Node Q, semi-major axis a, eccentricity e) were 
adjusted for optimum fit to the observations. The fitting pro- 
cedure utilizes a non-linear least-squares algorithm (Markwardt 
2005) that follows the direction of steepest descent of x 1 in the 
parameter space, x 1 being the reduced sum of squared devia- 
tions, i.e. the sum divided by the 13 degrees of freedom. To dis- 
tinguish between local and absolute minima, the initial parame- 
ters were varied over the broad ranges of their potential values: 
The semi major axis, a, was considered between 5.4 10 10 m and 
8.0 10 10 m, which corresponds to a total mass of Aa and Ab in 
the range 3 - 1OM . As specified in Sec.EJ e e [0.23,0.37]. 
The stellar diameters were examined between 0.4 and 12.4 mas. 
These limits refer respectively to the resolution limit of the in- 
terferometer and to the Roche lobe volume diameter D^. The 
latter is approximated to better than 1% by D^/d ~ 12.4mas 
(Eggleton 1983). If one of the stars were to have a diameter 
larger than Dl, the system would be an interacting binary and 
the simple model of two spherical, uniformly bright stars would 
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X =2.6 
a= (5.7 + -0.3)10 10 
D = 6.0 + -0.5 D 
D b = 3.3 + -0.6 D 
e= 0.230 + -0.005 
Q= 27.4 + -1.2 
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0.202 0.3565 



0.3570 
Orbital phase, 



0.377 



0.379 



0.937 



0.938 



Fig. 2. Corrected visibility values and standard deviation, confronted to a model of two uniformly luminous, spherical stars. The 
parameter values of the best fit (solid line) are indicated in the upper left corner. 



not apply. The position angle of the Ascending Node Q, mea- 
sured from North to East, equals if the projected orbital plane 
and the North-South axes are aligned. No previous measurement 
of Q. exists, and the angles D. and Q. + n can not be distinguished 
through interferometric measurements, Q. is therefore considered 
between and n. Varying the surface brightness ratio <p over 
the range specified in Eq.|8] has virtually no effect on y 2 , and 
<f> is therefore fixed at 1.46. Likewise, the period of the binary 
6 Vel(Aa-Ab) is fixed at T — 45. 150 days. No apsidal motion of 
the eclipsing system has been noted since its discovery in 2000. 
The orbital inclination has been fixed at 90°, although, given the 
stellar diameters and separations deduced in Section |2~2l the ac- 
tual inclination could lie between 87.5° and 92.5°. 

The best adjustment of the model to the measured visibil- 
ities and their 1-sigma statistical errors, is shown on Fig. [2] It 
corresponds to a reduced mean squared deviation x\ - 2.6, 
and is obtained for the following parameter values: a = (5.7 ± 
0.3) x 10 10 m,e = 0.230 ± 0.005, Q = (27.4 ± 1.2)°, D a = 
(6.0 ± 0.5) D @ ,D b = (3.3 ± 0.6) D e . The angle at primary 
eclipse is derived by the eccentricity as specified in Section lXTl 
a) = -(20 + 3)°. The parameter uncertainties equal the statistical 
errors, <x, scaled by the reduced mean deviation of the model to 
the measurements, i.e.^o cr. The dependence ofy 2 on the stellar 
diameters is illustrated on Fig. [3] 

The three visibilities measured on May 10, 2003 systemati- 
cally deviate from the model fit (see Fig.|2]i. There is no evident 
explanation for this deviation: the data were obtained under good 
atmospheric conditions and the calibrator was the same as on the 
other nights. If the three points are removed, the quality of the fit 
is improved, x\ = 1 -4, but within the uncertainties, the resulting 
parameter values are unchanged: a = (5.4 ± 0.5) x 10 I0 m, e = 
0.230 + 0.005, Q = (29.2 ± 2.4)°, D a = (6.6 ± 0.5) D ,D b = 
(3.2±O.5)£> . 

It is apparent from the relatively high x\\ that there are de- 
viations in addition to the purely statistical errors. They might 
be due to an underestimate of the calibrator's size or might re- 
flect some inaccuracies of the model of two uniformly bright, 
spherical stars. This is discussed in the subsequent section. 



4. Results and discussion 

4. 1 . The close eclipsing binary 6 Vel (Aa-Ab) 

The computations could be slightly biased if the diameter of the 
calibrator star were substantially misestimated, or if HD 63744 
were a - still undiscovered - binary system. On the other hand, 
HD 63744 is part of the catalog of interferometric calibrator stars 
by Borde et al. 2002, with its diameter (1 .63+0. 03)mas specified 
to a precision of 1.8%. Furthermore, it has been studied simul- 
taneously with other calibrator stars in VINCI observations by 
one of the authors (P. Kervella). In these investigations the vis- 
ibilities of HD 63744 equaled those expected for a single star of 
1 .63 ± 0.03 mas diameter. Thus, HD 63744 appears to be a reli- 
able calibrator. 

A perhaps more relevant aspect are possible astrophysical 
complexities of 5 Vel (Aa+ Ab) that are disregarded in the model 
of two uniformly bright, spherical stars. In particular, the rota- 
tional velocities of Aa and/or Ab are found to be high, with val- 
ues of -150 - 1 80km/ s (Royer et al. l2002l Hempel et al. [19981 
Holweger et al. 119991 ), which indicates that the two stars need 
not be uniformly luminous nor circular. 

Another possible over-simplification of our binary model 
is the constraint on the orbital inclination, i, being fixed at 
90°. Given the fitted semi-major axis and stellar diameters, we 
note that the eclipse durations (0.51 ± 0.05 days and 0.91 ± 
0.01 days) are shorter than they should be in the case of i = 90°, 
where the duration of the longer eclipse would have to exceed 
D a Tl{2na) = 1.06 days. We conclude that i is ~ 88° or ~ 92°, 
rather than 90°. All observations were performed out of eclipse 
and, therefore, the visibility values are nearly unaffected by such 
a small variation in i. With substantially more visibility measure- 
ments and an increased number of fitted parameters, the issue on 
the precise orbital inclination might be addressed in more detail. 

The most important and remarkable result of our analysis is 
that the stellar diameters of Aa and Ab are found to equal 6.0 ± 
0.5 D Q and 3.3 ± 0.6 D , respectively. This exceeds significantly, 
by factors ~ 1.4 - 3, the values expected if Aa and Ab are main 
sequence stars. If both diameters are constrained to lie below 
2.5D Q , the best fit corresponds to^ 2 = 16.7, which is far beyond 
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Tango, W. J., Davis, J., Thompson, R. J., Hanbury, R. 1979, PASAu, 3, 323 
Worley, C. E., & Douglass, G. G. 1997, A&AS, 125, 523 

Wright, C. O., Egan, M. P., Kraemer, K. E., Price, S. D. 2003 The Tycho-2 
Spectral Type Catalog 

Ever since the observations of Herschel (1847), cTVel has been 
taken to be a visual multiple star, <5Vel C and D being the outer 
components of the system. With my of 11.0 mag and 13.5 mag 
respectively (Jeffers et al. 1963), C and D would need to be 
of late spectral type, certainly no earlier than M, if they were 
at similar distance as <SVel (Aa+Ab+B). To our knowledge, the 
only existing spectra of C and D were recorded during a sur- 
vey of nearby M dwarfs (Hawley et al. 1996). While the limited 
range and resolution of the spectra precluded ready determina- 
tion of the spectral types of C and D, they were nevertheless 
estimated as ~G8V and ~K0V. Therefore, given their apparent 
magnitudes, C and D must be at much further distances than 
<5Vel (Aa+Ab+B). We conclude, that ciVel C and D are not phys- 
ically associated. Hence, (TVel ought to be classified as a triple 
stellar system only. 

5. Summary 

Seventeen VINCI visibility measurements of (5Vel(Aa+Ab) 
were fitted onto the model of two uniformly bright, spherical 
stars. The adjustment to the measurements does not provide in- 
dividual diameters compatible with A-type main sequence stars. 
The two stars appear, thus, to be in a more advanced evolution- 
ary stage. More data are, however, needed to ensure this result. 
As the stellar evolution is fast during this period, more detailed 
knowledge of the system might also constrain the models more 
tightly. Precise photometric and spectroscopic observations of 
the eclipses should provide the separate intensities and chemical 
compositions of Aa and Ab and, hence, permit further inferences 
on the age and evolutionary state of 5 Vel. 
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the present result and confirms that large diameters are required 
to account for the measured visibilities. 

4.2. On the physical association ofSVel C and D 



